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The regular structures obtained by optical lattice technology and their behaviour are analysed 
from the quantum information perspective. Initially, we demonstrate that a triangular optical lattice 
of two atomic species, bosonic or fermionic, can be employed to generate a variety of novel spin- 1/2 
models that include effective three-spin interactions. Such interactions can be employed to simulate 
specific one or two dimensional physical systems that are of particular interest for their condensed 
matter and entanglement properties. In particular, connections between the scaling behaviour of 
entanglement and the entanglement properties of closely spaced spins are drawn. Moreover, three- 
spin interactions are well suited to support quantum computing without the need to manipulate 
individual qubits. By employing Raman transitions or the interaction of the atomic electric dipole 
moment with magnetic field gradients, one can generate Hamiltonians that can be used for the 
physical implementation of geometrical or topological objects. This work serves as a review article 
that also includes many new results. 

PACS numbers: 03.75.Kk, 05.30.Jp, 42.50.-p, 73.43.-f 



I. INTRODUCTION 

With the development of optical lattice technology 
0- IE , considerable attention has been focused on the 
realisation of quantum computation 0, IE IE as weu 
as quantum simulation of a variety of many-particle sys- 
tems, such as spin chains and lattices [E IE El El This 
technology provides the possibility to probe and realise 
complex quantum models with unique properties in the 
laboratory. Examples that are of interest in various areas 
of physics are systems that include many-body interac- 
tions. These have been hard to study in the past due to 
the difficulty in controlling them externally and isolat- 
ing them from the environment |l2j . To overcome these 
proble ms, techniques have been developed in quantum 
optics pJlMlla which minimise imperfections and im- 
purities in the implementation of the desired structures, 
thus paving the way for the consideration of such "higher 
order" phenomena of multi-particle interactions. Their 
applications are of much interest to cold atom technol- 
ogy as well as to condensed matter physics and quantum 
information, some of which we shall see here. 

The initial point of our study is the presentation of the 
rich dynamics that governs the behaviour of an ultra cold 
atomic ensemble when it is superposed with appropriate 
optical lattices. For this purpose we consider the case 
of two species of atoms, denoted here by | and J. (see 
IE IE 111)')- trapped in the potential minima of a peri- 
odic lattice. These species can be two different hyperfine 
ground states of the same atom, coupled, via an excited 
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state, by a Raman transition. The system is brought 
initially into the Mott insulator phase where the num- 
ber of atoms at each site of the lattice is well defined. 
By restricting to the case of only one atom per site, it 
is possible to characterise the system by pseudo-spin ba- 
sis states provided by the internal ground states of the 
atom. Interactions between atoms in different sites are 
facilitated by virtual transitions. These are dictated by 
the tunneling coupling, J, from one site to its neigh- 
bours and by collisional couplings, U, that take place 
when two or more atoms are within the same site. Even- 
tually the evolution of the system can be effectively de- 
scribed by a wide set of spin interactions with coupling 
coefficients completely controlled by the tunneling and 
collisional couplings. This gives rise to the consideration 
of several applications that are mainly related with the 
three-spin interactions simulated on the lattice, which is 
the main focus of the present article. 

In that spirit, implementation of quantum simulation, 
of different physical models can be realised, with ground 
states that present a rich structure, such as multiple de- 
generacies and a variety of quantum phase transitions 
fltjl Ell Us) . Some of these multi-spin interactions have 
been theoreticallvstudied in the past in the context of the 
hard rod boson HEl2El2ll> usm g self-duality symmetries 
[22I. I2E]. Phase transitions between the corresponding 
ground states have been analysed 0> HH ■ Subsequently, 
these phases may also be viewed as possible phases of 
the initial system, in the Mott insulator, where the be- 
haviour of its ground state can be controlled at will [2E| • 
In this context the so-called cluster Hamiltonian is of 
considerable interest as its ground state exhibits unique 
entanglement properties. In this article we review the 
status quo of the existing analysis of these entanglement 
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properties as well as presenting new results that indicate 
interesting links between the entanglement properties of 
closely spaced spins and the scaling of entanglement with 
spin separation. 

To implement quantum computing, one can take ad- 
vantage of the three-spin interaction to construct multi- 
qubit gates that eventually can lead to quantum com- 
putation without the need to manipulate single qubits, 
referred to as 'global addressing'. We use a single qubit to 
localise operations, meaning that one and two qubit gates 
in a typical quantum computing scheme are replaced by 
two and three qubit gates in this scheme, which is what 
makes a triangular lattice with three-spin interactions 
such a natural environment for implementation of this 
concept. This global addressing lifts the stringent exper- 
imental requirement of single atom addressing for per- 
forming quantum gates. Moreover, error correction can 
be performed without the need to make measurements 
during the computation. 

The paper is organised as follows. In Section [HJ we 
present the physical system and the conditions required 
to obtain three-body interactions. The effective three- 
spin Hamiltonians for the case of bosonic or fermionic 
species of atoms in a system of three sites on a lattice 
are given in Sec. IIIII In Sec. IIVI we introduce com- 
plex couplings by considering the effect Raman transi- 
tions can have on the tunneling process and generalised 
effective Hamiltonians are presented that do not preserve 
the number of atoms of each species. In addition, the 
electric dipole of the atoms is considered which, through 
interaction with an external inhomogeneous field, can 
generate chirality. In Sec. IVIl a variety of entanglement 
properties of spin chains with three-spin interactions is 
analysed, and the cluster Hamiltonian is presented along 
with novel connections between the scaling behaviour of 
entanglement and the entanglement properties of closely 
spaced spins. A global addressing quantum computation 
model is presented in Sec. IVIII and, finally, in Sec. IVIIII 
concluding remarks are given. 

II. THE PHYSICAL MODEL 

Let us consider a cloud of ultra cold neutral atoms su- 
perimposed with several optical lattices 0, U 0, 0, |27j ■ 
For sufficiently strong intensities of the laser field, this 
system can be placed in the Mott insulator phase where 
the expectation value of only one particle per lattice 
site is energetically allowed |3|- This still allows for the 
implementation of non-trivial manipulations by virtual 
transitions that include energetically unfavourable states. 
Here, we are particularly interested in the setup of lat- 
tices that form an equilateral triangular configuration, as 
shown in Fig. ^ This allows for the simultaneous super- 
position of the positional wave functions of the atoms be- 
longing to the three sites. As we shall see in the following, 
this results in the generation of a three-spin interaction. 
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FIG. 1: The basic building block for the triangular lattice 
configuration. Three-spin interaction terms appear between 
sites 1, 2 and 3. For example, tunneling between 1 and 3 can 
happen through two different paths, directly and through site 
2. The latter results in an exchange interaction between 1 
and 3 that is influenced by the state of site 2. 

The main contributions to the dynamics of the atoms 
in the lattice sites are given by the collisions of the atoms 
within the same site and the tunneling transitions of the 
atoms between neighbouring sites. In particular, the cou- 
pling of the collisional interaction for atoms in the same 
site are taken to be very large in magnitude, while they 
are supposed to vanish when they are in different sites. 
Due to the low temperature of the system, this term is 
completely characterised by the s-wave scattering length. 
Furthermore, the overlap of the Wannier wave functions 
between adjacent sites determines the tunneling ampli- 
tude, J, of the atoms from one site to its neighbours. 
Here, the relative rate between the tunneling and the 
collisional interaction term is supposed to be very small, 
i.e. J <C U, so that the state of the system is mainly 
dominated by the collisional interaction. 

The Hamiltonian describing the three lattice sites with 
three atoms of species a — {Til} subject to the above 
interactions is given by 

H = H^ + V, (2.1) 

with 

Hi0) = \ u <y<y' a \a a \a> a i<y' a i<yi 

iaa' 

^ = -E( J >]^+i^+H.c.), 

where H.c. denotes the Hermitian conjugate, and aj a 
denotes the annihilation operator of atoms of species a at 
site j. The annihilation operator can describe fermions 
or bosons, satisfying commutation or anticommutation 
relations respectively, given by 

[aj<r,a i jl ,]± = 8j j >5 <T<T >, 

t t (2-2) 
[a ja ,af a >}± = [a) a ,a) l(r ,]± = 0, 

where the ± sign denotes the anticommutator or the com- 
mutator. The Hamiltonian is the lowest order in the 
expansion with respect to the tunneling interaction. 

Due to the large collisional couplings, activated when 
two or more atoms are present within the same site, the 
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weak tunneling transitions do not change the average 
number of atoms per site. This is achieved by adiabatic 
elimination of higher population states during the evolu- 
tion, leading to an effective Hamiltonian [2^ . The latter 
allows virtual transitions between these levels providing 
eventually non-trivial evolutions. According to this the 
low energy evolution of the bosonic or fermionic system, 
up to the third order in the tunneling interaction, is given 
by the effective Hamiltonian 
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The indices a, refer to states with one atom per site 
while 7, 8 refer to states with two or more atomic popu- 
lations per site, E 1 are the eigenvalues of the collisional 
part, H(°\ while we neglected fast rotating terms which 
is effective for long time intervals psj . 

It is instructive to estimate the energy scales involved 
in such a physical system. We would like to have a signif- 
icant effect of the three-spin interaction within the deco- 
herence times of the experimental system, which we can 
take here to be of the order of a few tens of ms. It is 
possible to vary the tunneling interactions from zero to 
some maximum value which we can take here to be of the 
order of J/fi -I kHz 0. In order to have a significant 
effect from the term J 3 /U 2 within the decoherence time, 
one should choose U/H ~10 kHz. This can be achieved 
experimentally by moving close to a Feshbach resonance 
[29L EfflL I3H ] , where U can take significantly large values. 
With respect to these parameters we have J/U ~ 10 _1 , 
which is within the Mott insulator regime, while the next 
order in perturbation theory is an order of magnitude 
smaller than the one considered here and hence negli- 
gible. This places the requirements of our proposal for 
detecting the effect of three-spin interactions within the 
range of the possible experimental values of the state of 
the art technology. 



III. THE EFFECTIVE THREE-SPIN 
INTERACTIONS 

The perturbative dynamics of the system is better pre- 
sented in terms of effective spin interactions. Indeed, 
within the regime of single atom occupancy per site, it 
is possible to switch to the pseudo-spin basis of states 
of the site j given by | f) = \rij-\ = = 0) and 

1 1) = = 0,rij[ = 1). Hence, the effective Hamilto- 
nian can be given in terms of Pauli matrices acting on 
states expressed in the pseudo-spin basis, as we shall see 
in the following. 



A. The bosonic model 

In this subsection and the following one, we will de- 
velop the Mott-Hubbard model in perturbation theory 



up to third order in J/U, i.e. the ratio of the tunnel- 
ing rate to the interaction term. In the first case, we 
will study a bosonic system, when two atoms of the same 
species are allowed to be in the same state. Eventually, 
our model is described by 



H cS = Y,[ A 3 l + B ^ 



3=1 
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A j 3 i+l + A 3 ^ 3 J+l + 3 + 



v (2) 



(3.1) 

where tr" is the a Pauli matrix at the site j. The cou- 
plings A, B, and \( % \ are given in terms of J a lU aa < by 
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where the symbol (t <_> 4) denotes the repeating of the 
same term as on its left, but with the j and J, indices 
interchanged. 

Knowing the dependence of the effective couplings al- 
lows us to modify the dynamics of the system at will, 
by changing the values of the tunneling rate or coupling 
constant, as seen in Fig. [5] Moreover, one body inter- 
actions in the Hamiltonian can be eliminated with an 
arbitrary Zeeman term of the form ■ B ■ cfj that can be 
added by applying a Raman transition with the appro- 
priate laser fields. One can also isolate different parts of 
Hamiltonian (|3.I(I , each one including a three-spin inter- 
action term, by varying the tunneling and/or the colli- 
sional couplings appropriately so that particular terms in 
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(|3.1(l vanish, while others are freely varied. An example 
of this can be seen in Fig. EJwhere the couplings A^ 1 ) and 
are depicted. There, for the special choice of the col- 
lisional terms, = Uu = 2.12Ufi, the A' 1 ) coupling is 
kept to zero for a wide range of the tunneling couplings, 
while the three-spin coupling, A*- 3 -* , can take any arbitrary 
value. One can also suppress the exchange interactions 
by keeping one of the two tunneling couplings zero, with- 
out affecting the freedom in obtaining arbitrary positive 
or negative values for A^ 3 \ as seen in Fig. 




FIG. 2: The effective couplings (a) A (1) , (b) A (2) , (c) A (3) and 
(d) A' 4 ) as functions of the tunneling couplings J T /[7 and 
J^/U. The tunneling couplings are set to be Jf = J% = Jf 
and the collisional couplings to be U-\-\ = Ufi = Uu = U. All 
the parameters are normalised with respect to U. 



B. The fermionic model 

In this subsection, we will show the effective Hamilto- 
nian for a system where two atoms of the same species 
are not allowed to be in the same site, so that they are 
described by fermionic operators. For the same reason, 
there is just one collisional coupling U-\± = U, and the 
others can be thought to contribute with an infinite en- 
ergy, U^,Un — > oo. After a tedious calculation and 
keeping terms to third order in Jf /U, the effective Hamil- 
tonian appears as, 



H, 



off : 



<T 3 ,J J + 



!)+M (3 V 



'1 U 2 U 3/ 
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FIG. 3: The effective couplings A (1) and A (3) are plotted 
against J T /U and J 1 /U for U n = U u = 2.12C/ and U u = U. 
The coupling A' 1 ' appears almost constant and zero as the 
unequal collisional terms can create a plateau area for a small 
range of the tunneling couplings, while A' 3 - 1 can be varied 
freely to positive or negative values. 



where the effective couplings are functions of the tunnel- 
ing and collisional couplings, given by 
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In this case, the dependence of the coupling terms on 
the parameters of the initial Hamiltonian is simpler than 
in the bosonic one. If the tunneling constants do not 
depend on the pseudo-spin orientation, indicated here by 
the subscript j = 1,2,3, then any three-spin interaction 
vanishes. Nevertheless, when the tunneling amplitudes 
depend on the spin and there is just one of the orientation 
with non-zero tunneling, then, only the two- and three- 
spin interactions in the z direction remain. A general 
picture of their behaviour can been seen in Fig. 21 



IV. COMPLEX COUPLINGS 
A. Raman activated tunneling 

From the previous models, it is possible to create new, 
different, Hamiltonians by employi ng techniques avail- 
able from quantum optics H, 0, |32] , like the application 
of Raman transitions during the tunneling process. If 
the lasers producing the Raman transition form stand- 
ing waves, it is possible to activate tunneling transitions 
of atoms that simultaneously experience a change in their 
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FIG. 4: The effective couplings (a) (b) fi (2) , (c) ^ (3) and 
(d) M < 4 > as functions of the tunneling couplings J T /U and 
J^/U, where the tunneling couplings are set to be Jf = J% = 

J 3 ■ 

internal state. As we shall see in the following, the result- 
ing Hamiltonian is given by an SU(2) rotation applied to 
each Pauli matrix of the previous Hamiltonians. 

Consider the case of activating the tunneling with the 
application of two individual Raman transitions. These 
transitions consist of four paired laser beams Li, L2 and 
L[, L' 2l each pair having a blue de-tuning A and A', 
different for the two different transitions. The phases 
and amplitudes of the laser beams can be properly tuned 
so that the Raman transitions allow the tunneling of the 
states 

|+) =cos6»|a) +sin0e- i *|&) 
|-) = sin0|a} -cos0 e - # |&), 
or in a compact notation, 

(II)-**"®) <«> 

with the unitary SU (2) matrix 

/ 1 n\ ( cos 6 e 1 ^ sin 6 \ , . 

^={sin9 -e*>co S 6) ^ 

In the above equations, 4> denotes the phase difference 
between the Li laser field, while tan 9 = l^/^i |- 
are the Rabi frequencies of the laser fields. Hence, the 
resulting tunneling Hamiltonian can be obtained from 
the initial one via an SU(2) rotation, 

Vc = 9 Vg^ = - J2(J+cf4 +1 + J-C^ c r +1 + H.c). 

i 

where the corresponding tunneling couplings are formally 
identified, i.e. J + = ,P and J~ = J*. Note that the col- 
lisional Hamiltonian is not affected by the Raman transi- 
tions, and hence it is not transformed under g rotations. 




FIG. 5: Example of a Raman transition activated by a pair 
of blue detuning laser fields L\ and L2. 



It is easy to derive the effective Hamiltonian for this 
transformation using the perturbative expansion. From 
(|2.3() we straightforwardly obtain the second, third and, 
likewise, any order term of the Hamiltonian H c ^ that ap- 
pear after the application of the Raman transition. They 
are given by an SU(2) rotation that acts on the Pauli 
matrices of the initial effective Hamiltonian and read 

H$(cb,e)=g(cb,e)H^gHcb,e), 

where n is the order of the perturbation. 

This approach provides a variety of control parameters 
(e.g. the angle <f> an d the ratio of the couplings of the two 
added Hamiltonians) and, in addition, one can have these 
variables independent for each of the three directions of 
the two dimensional (triangular) optical lattice. Particu- 
lar settings of these structures have been proved to gen- 
erate topological phenomena that support exotic an- 
ionic excitations, useful for the construction of topologi- 
cal memories [33] ■ In addition, the possibility of varying, 
arbitrarily, the control parameters gives us the natural 
setup to study such phenomena as geometrical phases in 
lattice systems. 

B. Complex tunneling and topological effects 

Consider the case where we employ complex tunnel- 
ing couplings |3j[ to the optical lattice evolution. This 
can be performed by employing additional characteris- 
tics of the atoms, like an electric moment d e and an 
external electromagnetic field. As the external field 
can break time reversal symmetry, new terms of the 
form {crj crJ +1 o-j : +2 — o']'ctJ +1 ct| +2 } appear in the effec- 
tive Hamiltonian. In particular, the minimal coupling 
deduced from the electric dipole of the atom with the 
external field can be given, in general, by substituting its 
momentum by 

p^p + (d e - V)A(x), 

where A is the corresponding vector potential. The new 
term satisfies the Gauss gauge if we demand that V • A = 
0, hence it can generate a possible phase factor for the 
tunneling couplings. 
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Due to an evolution dictated by the differential form 
of the Aharonov-Bohm effect j3f| the cyclic move of an 
electric moment through a gradient of a magnetic field 
contributes the phase 



(d e ■ V)A ■ ds, 



(4.3) 



to the initial state, where S is the surface enclosed by 
the cyclic path of the electric moment. By inspection 
of relation (|4.3|) we see that a nontrivial phase can be 
produced if we generate an inhomogeneous magnetic field 
in the neighbourhood of the dipole. In particular, a non- 
zero gradient of the magnetic component perpendicular 
to the surface S, varying in the direction of the dipole, 
ensures a non-zero phase factor. For example, if we take 
S to lie on the x-y plane and d e is perpendicular to the 
surface 5, then a non-zero phase, </>, is produced if there 
is a non-vanishing gradient of the magnetic field along 
the z direction. This is sketched in Figure HHa), where 
the magnetic lines are plotted such that they produce the 
proper variation of B z in the z direction. Alternatively, 
if d e is along the surface plane, then a non-zero phase is 
produced if the z component of the magnetic field has a 
non- vanishing gradient along the direction of d e as seen in 
FigureEJb), where only the z component of the magnetic 
field has been depicted. 
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FIG. 6: The path circulation of the electric dipole in the 
inhomogeneous magnetic field. Figure (a) depicts magnetic 
field lines sufficient to produce the appropriate non- vanishing 
gradient of B z along the z axis. Figure (b) depicts only B z 
and how it varies along the direction of the dipole, d e . 

We can denote by J = e l ^\ J\ the generation of the 
phase factor of the tunneling coupling between two neigh- 
bouring sites, with 



Xi + l 



(d e ■ V)A ■ dx. 



Here Xi and Xi+\ denote the positions of the lattice sites 
connected by the tunneling coupling J. 

In order to isolate the new effects generated by the 
consideration of complex tunneling couplings, we restrict 
ourselves to purely imaginary ones, i.e. Jf = ±i\ J?\, We 
also focus, initially, on the case where the optical lattices 
generate a two dimensional structure of equilateral trian- 
gles, as in Figure[7| Such a non-bipartite structure is nec- 
essary in order to manifest the breaking of the symmetry 
under time reversal, T, in our model, eventually produc- 
ing an effective Hamiltonian that is not invariant under 



complex conjugation of the tunneling couplings. More- 
over, as the second order perturbation theory is mani- 
festly T symmetric, we need to consider the third order. 
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FIG. 7: The basic building block for the triangular lattice 
configuration. Exchange of atoms through sites i and j can 
happen directly or through site k, but with a phase difference, 
causing the generation of the chiral three-spin interaction. 

In this case, the effective Hamiltonian |3y|, up to the 
second order in perturbation theory, becomes 



rr(2) _ 
H eff ~ 



n 



Al + Baf 



+ D(a*a* +1 +aia» 



i+i. 



The above couplings are given in perturbation theory by 
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It is easily seen that the values of the above interac- 
tions remain unchanged compared to the case with no 
magnetic field. This is the case since, up to the sec- 
ond order perturbation theory, there are no contributions 
from circular paths that can experience the magnetic 
field. Note that it is possible to apply Raman transi- 
tions that exactly cancel the a z Zeeman term. Also, one 
can choose U-\i — > oo (possibly by Feshbach resonances) 
so that C vanishes. In addition, one of the collisional 
couplings can be chosen to be attractive, satisfying e.g. 
f/pr = —Uu = —U, so that, for J T = , the A cou- 
pling vanishes. Such attractive collisional couplings can 
be achieved by Feshbach resonances [32J • Stability of the 
system is maintained if the attractive couplings are acti- 
vated adiabatically, while the atomic ensemble is in the 
Mott insulator regime. Since attractive couplings may 
eventually lead to pairing on the same lattice site, this 
negative U interaction can only be applied for a short 
period. 

The third order in the perturbative expansion includes 
interactions between neighbouring sites via the circular 
path around the elementary triangle. This brings in the 
effect of the magnetic field, contributing a phase factor. 
While new interaction terms appear, the couplings A, 
B, C and D remain the same. Considering the case of 
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collisional and tunneling couplings, such that C = D = 0, 
the effective Hamiltonian becomes, up to global Zeeman 
terms, 



H. 



(3) 



£ [E{p?J> - o\o*) + Fe lmn a\afal] 
(ijk) 



where 



7T jl 



/T jl 

J 1 ), F = z^(jT_jl). 



Here (zj'fc) denotes nearest neighbours, ei mn with 
{I, m, n} = {#, y, z} denotes the total antisymmetric ten- 
sor in three dimensions and summation over the indices 
l,m,n is implied. From the expression of the effective 

(3) 

Hamiltonian, H^J, it is apparent that the E and F cou- 
plings are not invariant under time reversal, T, that is, 
they change sign after complex conjugation. This leads 
to the breaking of the chiral symmetry between the two 
opposite circulations that atoms can take around a tri- 
angle. 

By additionally taking = — = J, one can set all 
the couplings to be zero apart from F and the effective 
Hamiltonian reduces to 



(ijk) 



(4.4) 



with a = {a x ,a v ,a z ) and F = \J\ 3 /U 2 . Remarkably, 
with this physical proposal, the interaction term (|4.4|) 
can be isolated, especially from the Zeeman terms that 
are predominant in equivalent solid state implementa- 
tions. This interaction term is also known in the liter- 
ature as the chirality operator |38| . It breaks the time 
reversal symmetry of the system, a consequence of the 
externally applied field, by effectively splitting the de- 
generacy of the ground state into two orthogonal sectors, 
namely "+" and "— ", related by time reversal, T. These 

sectors are uniquely described by the eigenstates of H cS ' 
at the sites of one triangle. The lowest energy sector with 
eigenenergy E+ = ~2y/3F is given by 



1 



l*i/ 2 ) = 7^0 TU> +"\ TIT) +^ 2 UTT» 



I* 



-1/2/ 



V3 



(1 ut>+H m>+^ 2 i m» 



(4.5) 



where 



,3 - 



= 1. The excited sector, 1^ 



±1/2/ 



repre- 
sents counter propagation with eigenvalue £L — 2\^3F 
and it is obtained from (|4.5|l by complex conjugation 
HE HE El • It has been argued that this configuration, a 
result of frustration due to the triangular lattice, and the 
disorder due to the presence of the magnetic field, leads 
to analogous behaviour with the fractional quantum Hall 
effect HE EJ and, in particular, it can be described by 
the m = 2 Laughlin wavefunction defined on the lattice 
sites. A recent example demonstrating this analogy is 
given in |4^ . 



V. ONE- AND TWO-DIMENSIONAL MODELS 

It is also possible to employ the three-spin interactions 
that we studied in the previous sections for the construc- 
tion of extended one and two dimensional systems. The 
two dimensional generalisation is rather straightforward 
as the triangular system we considered is already defined 
on the plane. Hence, all the interactions considered so 
far can be generalised for the case of a two dimensional 
lattice where the summation runs through all the lattice 
sites with each site having six neighbours. 

The construction of the one dimensional model is more 
involved. In particular, we now consider a whole chain 
of triangles in the one dimensional pattern shown in Fig. 
[H] In principle this configuration can extend our model 
from the triangle to a chain. Nevertheless, a careful con- 
sideration of the two spin interactions shows that terms 
of the form o~faf, 2 appear in the effective Hamiltonian, 
due to the triangular setting (see Fig. [SJ. Such Hamilto- 
nian terms involving nearest and next-to-nearest neigh- 
bour interactions are of interest in their own right p^lllSf 
but will not be addressed here. It is possible to introduce 
a longitudinal optical lattice with half of the initial wave- 
length, and an appropriate amplitude such that it cancels 
exactly those interactions generating, finally, chains with 
only neighbouring couplings (for more details, see lVll D"|l . 




1 



FIG. 8: The one dimensional chain constructed out of equi- 
lateral triangles. Each triangle experiences the three-spin in- 
teractions presented in the previous sections. 



In a similar fashion it is possible to avoid generation 
of terms of the form crfaf +2 + of of +2 by deactivating 
the longitudinal tunneling coupling in one of the modes, 
e.g. the | mode, which deactivates the corresponding 
exchange interaction. 

As we are particularly interested in three-spin in- 
teractions we would like to isolate the chain term 

Ei( cr f<+ i cr f+2 + °f Qj-tW+a)' the term of Hamilto- 
nian This term includes, in addition, all the pos- 
sible triangular permutations. To achieve that, we could 
deactivate the non-longitudinal tunneling for one of the 
two modes, e.g. the one that traps the "f atoms. The in- 
teraction of crf +1 of +2 is homogeneous, hence it does not 
pose such a problem when it is extended to the one di- 
mensional ladder. With the above procedures, we can 
finally obtain a chain Hamiltonian as in where the 
summation runs up to the total number, N, of lattice 
sites. 
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This class of Hamiltonians gives rise to a rich variety 
of ground states for the spin ladder. Let us consider, as 
an example, the case when the tunnelling for one of the 
species is reduced to zero. This switches off all terms 
involving <j x and cr y in the Hamiltonian (|3 -II) . Suppose 
further that we have used additional laser fields to can- 
cel out the Zeeman term Ba z . Then, the only non-zero 
coefficients in the Hamiltonian (|3.1|) are X^\ the Ising 
interaction, and A' 3 -* , the 3-spin interaction. The Hamil- 
tonian becomes 



H = \ ( > 1 



'i+l 



(5.1) 



Without the 3-spin interaction, the Hamiltonian de- 
scribes the classical Ising model. The sign of the Ising 
interaction can be tuned by the relative magnitudes of the 
repulsive potentials U. For negative \^\ the system is 
ferromagnetic. The up-down symmetry in the z-direction 
of the spin can be spontaneously broken and there are two 
degenerate ground states: all spins up, or all spins down. 
For positive A' 1 ), the ground state is the antiferromag- 
netic Neel state, | UUU ■■■) or | UUU ■■■)■ Note 
that this breaks, spontaneously, the lattice translational 
symmetry with a spin configuration of lattice period 2. 

On the other hand, if we switched off the Ising interac- 
tion with a suitable choice of t/'s, the 3-spin term would 
give rise to different broken-symmetry ground states. 
For A^ 3 ' < 0, there are four degenerate classical ground 
states: a ferromagnetic state | TU ■ * * ) and three states 
with lattice period 3, | TUTU ■■•>, I ITUT1 ■■■) and 
I UTUT •••)■ For ^ (3) > 0) the ground states are the 
same with a z reversed at each site. Note that the 3- 
spin term explicitly breaks spin-reversal symmetry, re- 
flecting the difference in the hopping amplitudes of the 
two atomic species of the original system. Nevertheless, 
the lattice translational symmetry is spontaneously bro- 
ken in the period-3 states. A small Zeeman field in the 
^-direction or a small antiferromagnetic Ising coupling 
will select out the period-3 states from the ferromagnetic 
state. 

We can switch directly between the period-2 and 
period-3 ground states by tuning the interaction parame- 
ters A^ 1 ) and X^ . The period-2 Neel state has an energy 
of — A^ 1 ' per site while the lowest-energy period-3 states 
have an energy of — AW /3—|A^I per site. So, we expect a 
first-order phase transition along the line 2A^ = 3|A^ 3 ^|. 

We can introduce quantum correlations into these 
models by introducing a transverse Zeeman field, B x a x , 
in the Hamiltonian at each site 



H = 



B x ai 



(5.2) 



This transverse field gives rise to nipping between the 
up and down spin states (of the a z basis). The Ising 
chain (non-zero X^ with A*- 3 -* = 0) in a transverse field 
is a well-known model. This is discussed extensively in 
(T?| and it can be analytically solved using the Jordan- 
Wigner transformation. As we increase the transverse 



field B x from zero, the magnetisation in the z-direction 
(ferromagnetic or antiferromagnetic according to the sign 
of A^ 1 )) is reduced. It vanishes at the quantum critical 
point B x = \X^ 1 '\. At larger values of the transverse 
field, the spins develops a polarisation in the a;-direction. 
This quantum phase transition for the one-dimensional 
chain belongs to the same universality class as the two- 
dimensional classical Ising model with no transverse field. 

What happens when we introduce a weak 3-spin inter- 
action into the Ising chain with transverse field? Unfor- 
tunately, the model is not exactly solvable. As already 
mentioned, the 3-spin term breaks spin-reversal symme- 
try and so, just like a Zeeman field in the z-direction, it 
is a relevant perturbation in the renormalisation group 
sense. 

Indeed, considering the opposite limit where the 3-spin 
interaction is strong, we know that a phase transition ex- 
ists in the presence of a transverse field with the 3-spin 
interaction but no 2-spin Ising interaction |22J, |23j . At 
this transition, the period-3 order parameter vanishes at 
B x = lA^I- Although this transition looks analogous to 
the one found in an Ising chain in a transverse field, it 
does not belong the same universality class as the Ising 
chain. It is believed that it belongs to the same univer- 
sality class as the classical two-dimensional 4-state Potts 
model [H|23. 

This leads us to speculate about how the direct tran- 
sition between period-2 and period-3 ground states is 
affected by the quantum fluctuations introduced by a 
transverse field. From the viewpoint of condensed mat- 
ter physics, it will be interesting to study the critical 
behaviour of this model. We believe that the critical 
point belongs to the same universality class as the two- 
dimensional 3-state Potts model, borrowing arguments 
described in |43| for the melting of commensurate struc- 
tures. 

Furthermore, we note that the phases of this spin chain 
are similar to the quantum hard-rod system studied in 
[l8|. That model has the additional complexity of a 
macroscopically large number of classically degenerate 
ground states at the classical transition (without trans- 
verse field) . We believe that our system is easier to imple- 
ment in the context of optical lattices for neutral atoms. 

From the viewpoint of quantum information, it is in- 
teresting to study the degree of entanglement in this sys- 
tem. Indeed, we can define a reduced density matrix for 
L contiguous spins in a system with N spins by tracing 
out the other spins of the system, 

PL = trjv-il^oX^ol 

where \^o) is the ground state of the system (4^. A 
measure of how these L spins are entangled with the rest 
of the chain is the von Neumann entropy, 

Sl = -tr(p L log 2 Pl) 

We find that the entanglement between two halves of the 
system increases dramatically near the critical point (Fig. 
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|5J) , similar to other studies of systems near criticality . 
It will be interesting to investigate whether the period- 
2 and period-3 states can be entangled by driving the 
system through the transition. 




FIG. 9: The entropy Sl as a function of the length L of the 
reduced density matrix and the transverse magnetic field B x 
for A' 3 ' = 1. The critical behaviour is indicated at the point 
B* « 1. 



VI. LOCALISABLE ENTANGLEMENT AND 
CLUSTER HAMILTONIAN 

In the past, Hamiltonians describing thre e-sp in inter- 
actions have been of limited interest Hj as they 
were difficult to implement and control experimentally. 
The previous results demonstrate that Hamiltonians with 
three-spin interactions can be implemented and parame- 
ters controlled across a wide intervals. One may suspect 
that ground states of three-spin interaction Hamiltonians 
exhibit unique properties as compared to ground states 
generated merely by two-spin interactions. This moti- 
vates the study of the properties of the ground state of a 
particular three-spin Hamiltonian for different paramet- 
ric regimes. Possible phase transitions induced by vary- 
ing these parameters are explored employing two possible 
signatures of critical behaviour that are quite different in 
nature. In particular, new critical phenomena in three- 
spin Hamiltonians that cannot be detected on the level 
of classical correlations will be demonstrated. 

(i) A traditional approach to criticality of the ground 
state studies two-point correlation functions between 
spins 1 and L, given by Cf/ = (afof) - (of) (of), for 
varying L, where a, j3 — x,y,z. These two-point correla- 
tions may exhibit two types of generic behaviour, namely 
(a) exponential decay in L, i.e. the correlation length £, 
defined as 



r 



i 



lim y log Cll , 



(6.1) 



is finite or, (b), power-law decay in L, i.e. ~ L~ q 

for some q, which implies an infinite correlation length £ 
indicating a critical point in the system . 

(ii) While the two-point correlation functions C^f* are 
a possible indicator for critical behaviour, they provide 
an incomplete view of the quantum correlations between 
spins 1 and L. Indeed, they ignore correlations through 
all the other spins, by tracing them out. Considering, for 
example, the GHZ state, \GHZ) = (|000) + |H1))/V% al- 
ready shows that this loses important information. Trac- 
ing out particle 2 leaves particles 1 and 3 in an unen- 
tangled state. However, measuring the second particle in 
the cr^-eigenbasis leaves particles 1 and 3 in a maximally 
entangled state. Therefore one may define the localisable 
entanglement between spins 1 and L as the largest 

average entanglement that can be obtained by perform- 
ing optimised local measurements on all the other spins 
|45|. In analogy to Eqn. 
glement length, £ E , by 

G 1 



H6.1fl one can define the entan- 



lim — log EY 



(loc) 



(6.2) 



It is now an interesting question whether criticality ac- 
cording to one of these indicators implies criticality ac- 
cording to the other. The localisable entanglement length 
is always larger than, or equal to, the two-point correla- 
tion length and indeed, it has been shown that there are 
cases where criticality behaviour can be revealed only by 
the diverging localisable entanglement length while the 
classical correlation length remains finite [lfl l4fH . Such 
behaviour is also expected to appear when we consider 
particular three-spin interaction Hamiltonians. To see 
this consider the Hamiltonian 



H = 



'i+l 



Baf 



(6.3) 



where we assume periodic boundary conditions. The fact 
that <jf_i<J*: af +1 commute for different i and employing 
raising operator LI = a 



%a 



allows to deter- 



mine the entire spectrum of H for B = quite easily. 
The unique ground state of H for B = is the well- 
known cluster state |4tI ligj , which has previously been 
studied as a resource in the context of quantum compu- 
tation. It possesses a finite energy gap of AE — 2 above 
its ground state jl^ ■ For finite B the energy eigenvalues 
of the system can still be found using the Jordan- Wigner 
transformation and a lengthy, but straightforward, cal- 
culation shows that the energy gap persists for \B\ ^ 1. 
The exact solution also shows that the system has criti- 
cal points for \B\ = 1 at which the two-point correlation 
length and the entanglement length diverge. For any 
other value of B and in particular for B = 0, the sys- 
tem does not exhibit a diverging two-point correlation 
length as is expected from the finite energy gap above 
the ground state. Indeed, correlation functions such as 



(6.4) 
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with 



i>ab 



and 



1 



4tt J_ 27r ^/B 2 + 1 + 2B cos r 



(b-a)r 
■. sm ar, 



Xab 



1 

47T 



(b — a)r 



-2k \AB 2 + 1 + 2B cos r 



(6.5) 



dr. 



(6.6) 

can be computed and the corresponding correlation 
length can be determined analytically using standard 
techniques (see e.g. Fig. (fTU|l ) |50|. The two-point cor- 
relation functions, such as Eqn. I|6.4(l . exhibit a power- 
law decay at the critical points \B\ = 1 while they de- 
cay exponentially for all other values of B in contrast to 
the anisotropic XY-model whose Cff correlation func- 
tion tends to a finite constant in the limit of L — ► oo 
for \B\ < 1 [53. This discrepancy is due to the finite 
energy gap the model in Eqn. I|6.3|) exhibits above a 
non-degenerate ground state in the interval \B\ < 1. 

When we study three-spin interactions it is natural 
to consider the behaviour of higher-order correlations. 
For the ground state with magnetic field B = 0, all 
three-point correlation except, obviously, (of—icf af + i) 
vanish. Indeed, if we consider n > 4 neighbouring 
sites and choose for each of these randomly one of the 
operators u x ,a y ,cr z or I then the probability that the 
resulting correlation will be non-vanishing is given by 
p = 2~( 2+ "). For \B\ > 0, however, far more corre- 
lations are non-vanishing and the rate of non-vanishing 
correlations scales approximately as 0.858™. This marked 
difference, which distinguishes B = 0, is due to the higher 
symmetry that the Hamiltonian exhibits at that point. 

In the following we shall consider the localisable en- 
tanglement and the corresponding length as described in 
(ii). Compared to the two-point correlations, the com- 
putation of the localisable entanglement is considerably 
more involved due to the optimisation process. Never- 
theless, it is easy to show that the entanglement length 
diverges for B = 0. In that case the ground state of the 
Hamiltonian (|6.3|) is a cluster state with the property 
that any two spins can be made deterministically max- 
imally entangled by measuring the a z operator on each 
spin in between the target spins, while measuring the o~ x 
operator on the remaining spins. Indeed, this property 
underlies its importance for quantum computation as it 
allows us to propagate a quantum computation through 
the lattice via local measurements |4gj . 

For finite values of B it is difficult to obtain the exact 
value of the localisable entanglement. Nevertheless, to 
establish a diverging entanglement length it is sufficient 
to provide lower bounds that can be obtained by pre- 
scribing specific measurement schemes. Indeed, for the 
ground state of i|6.3[l in the interval \B\ < 1, consider 
two spins 1 and L = 2k + 1 where k € N. Measure 
the o~ x operator on spin 2 and on all remaining spins, 
other than 1 and L, measure the o~ z operator. By know- 
ing the analytic form of the ground state one can ob- 



tain the average entanglement over all possible measure- 
ment outcomes in terms of the concurrence, that tends 

to E'oc = (l — l-Bl 2 ) 1 ^ 4 for k — ► 00. This demonstrates 
that the localisable entanglement length is infinite in the 
full interval \B\ < 1. This surprising critical behaviour 
for the whole interval \B\ < 1 is not evident from simple 
two-point or n-point correlation function which exhibit 
finite correlation lengths. For \B\ > 1 however, numer- 
ical results, employing a simulated annealing technique 
to find the optimal measurement for a chain of 16 spins, 
show that the localisable entanglement exhibits a finite 
length scale. 



Entanglement length 

- - Correlation length 





1 Magnetic Field B 1 



FIG. 10: Both, the two-point correlation length for C" 
(dashed line) and the localisable entanglement length (solid 
line) are shown for various magnetic field for chain of length 
16. Note that the localisable entanglement length diverges 
in the whole interval \B\ < 1 while the two-point correlation 
length is finite in this interval. 



From a practical point of view it is of interest to see 
how resilient these predicted features are with respect 
to noise. To this end, note that in Fig. I|10|) both the 
two-point correlation length and localisable entanglement 
length are drawn versus the magnetic field. In the in- 
terval \B\ < 1 the entanglement length diverges while 
the correlation length remains finite. For finite tempera- 
tures the localisable entanglement becomes finite every- 
where but, for temperatures that are much smaller than 
the gap above the ground state, it remains considerably 
larger than the classical correlation length. This demon- 
strates the resilience of this phenomenon against thermal 
perturbations. 

The preceding analysis of our three-spin interaction 
Hamiltonian 16. 3|) has indicated a marked qualitative 
discrepancy between the localisable entanglement length 
and the two-point correlation length with regard to varia- 
tions in the magnetic field strength parameter, B: across 
the interval between the two critical points, \B\ < 1, 
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the former quantity diverges while the latter remains fi- 
nite. We surmised that knowledge of simple two-point 
correlation functions (such as C*£) alone is insufficient 
in order to arrive at a complete characterisation of crit- 
icality. Nevertheless, the question remains whether the 
aforementioned discrepancy could not be reconciled by 
some non-trivial combination (motivated by quantum in- 
formation theory) of two-point correlation functions. A 
natural and easily computable candidate for such a quan- 
tity is the logarithmic negativity [&l| (for basic ideas of 
the theory of entanglement measures see, for example, 
|52|). which serves as a practical tool in the endeavour 
to quantify multipartite entanglement. The logarithmic 
negativity of a quantum state p is defined as 

E N = log 2 tr|p r |, (6.7) 

where p T denotes the partial transpose of the density 
operator and tr\.\ refers to the trace- norm, i.e. the sum 
of the singular values of the operator. 

We are interested in the entanglement, as measured by 
the logarithmic negativity, between two spin-i particles 
residing at sites i and j on the chain modelled by Hamil- 
tonian l|(j.3[l . The composite state for the pair of spins 
is described by the two-site reduced density matrix, p^, 
which is formally obtained from the total density matrix 
(describing the ground state of the whole chain) by trac- 
ing out all the spin degrees of freedom apart from the 
two spins under consideration. 

In practice, however, the easiest procedure is to ex- 
pand the general two-site density matrix in the trace- 
orthogonal basis formed by the tensor products of the 
Pauli spin operators at either site as 

a, 0=0 

where the expansion coefficients are expectation values 
with respect to the ground state of the Hamiltonian. It 
is useful to represent this operator expansion in the stan- 
dard two-qubit basis, {|TT) , |tl) , lit) , III)}, where |T> 
and ||) represent the eigenstates of o~ z with eigenvalues 
+1 and —1, respectively. In that basis, 

a,P=0 

withs = I TXT ] = + .s 1 = 10(11 = 
s 2 ^IT>UI = <7+ s 3 ^ll}(TI = ^. 

A priori, the construction of pij then requires knowl- 
edge of sixteen expansion coefficients, (sfsj). Owing to 
symmetry properties of the Hamiltonian, however, these 
are not all independent. For example, the Hamiltonian 
possesses the global phase- flip symmetry [H,U] = 0, 
where U = ®f=\0~ z . This symmetry then carries over 
directly to the two-site reduced density matrix obtained 
from the global ground state in the form [pij, cfcfj] = 0, 



forcing half of its matrix elements to vanish. Further- 
more, translational invariance of the Hamiltonian dic- 
tates that Pi i,n = Pii.ii (here and henceforth wc drop 
the indices i and j). Finally, the reduced density matrix 
must be a real, symmetric matrix with unity trace. 

Having taken the Hamiltonian's symmetries into con- 
sideration, we find that the reduced density matrix pos- 
sesses three distinct diagonal elements, 

Pi = PTT.TT = (1 + 2(^~) + (o?oJ»/4, 

P2 = P1U11 = C 1 - (°f°j))/ 4 : 

Ps = PlU1 = (l-2(a z ) + { ( T?ap)/4, 
and two distinct non-zero off-diagonal elements, 

P+ = PTT.U = P- = PTUT = ( a t a j)- 

The reduced density matrix is therefore uniquely deter- 
mined by the set of four expectation values, {(a z ) ,(cr| <rj) , 

Since the reduced density matrix between nearest 
neighbours does not exhibit entanglement for any value 
of the magnetic field B we will now consider the reduced 
density matrix with respect to a "bridge pair" of spins 
(any pair of spins sandwiching a third) . Given the trans- 
lational invariance of the Hamiltonian, we may restrict 
our study to sites 1 and 3, without loss of generality. 
The partial transpose of the reduced density matrix with 
respect to cither subsystem differs from the original ma- 
trix only through an interchange of p+ and p_ (along 
with {i,j} and {1,3}). Its four eigenvalues are given by 



(pi + p 3 ) ± (pi + p 3 ) 2 - 4(pip 3 - P-) 



and 

^3,4 = 92 ± P+- 

In order to quantify the entanglement within our 
bridge pair, all that remains is to compute the set of 
four simple expectation values that constitute the raw 
ingredients for the pair's reduced density matrix. We are 
especially interested in the thermodynamic limit, where 
N — > oo and we are no longer dealing with complicated 
sums, but with manageable integrals. In fact, all four 
expectation values are simple combinations of the two 
integrals, Eqs. iltj.5fl and i|tj.(jfl . defined earlier: 

(c rZ )=-Xoo, Ko-f)=Xoo+^i3-Xi3, 

K4> - ^13 (O, ) = lxi3(<? Z ) ■ 

For a few special values of B, these integrals, and the 
associated logarithmic negativity, can be computed pre- 
cisely with ease. The results are summarised in the table 
below. 
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logarithmic negativity 



For general B, the integrals are most easily computed 
numerically. A plot of the logarithmic negativity versus 
the magnetic field strength is shown in Fig. i|ll|) . 



0.3 




o 



0.05 - 

o- ' 1 - 

-5 -4 -3 -2 -1 1 2 3 4 5 

Magnetic Field B 

FIG. 11: The logarithmic negativity is shown as a function 
of the magnetic field strength parameter B in the thermody- 
namic limit TV — > oo. 



From among the features discernible in Fig. Qllfl. we 
are primarily interested in the interval enclosed by the 
two critical points, \B\ ^ 1, where the logarithmic nega- 
tivity appears to vanish identically. That characteristic 
would be directly analogous to the diverging localisable 
entanglement length encountered in Fig. 1)10(1. Close in- 
spection reveals, however, that the logarithmic negativ- 
ity grows positive even before reaching the critical points, 
B c = ±1. This phenomenon already suggested itself as 
a result of our analytical computation of the logarith- 
mic negativity (c.f. the provided table), and is most dis- 
cernibly manifest in Fig. i|12|) . which provides a close-up 
of the area immediately surrounding B c = 1. 
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FIG. 12: Here, Fig. HI IB is magnified in a small interval 
centred around the point B c — 1. The logarithmic negativity 
is clearly seen to take on positive values before reaching B c = 
1. 

Note that we may bypass many steps involved in the 
computation of the logarithmic negativity if we are not 
concerned with its precise value across the whole range 
of B, but are merely interested in establishing the region 
where it vanishes. For that purpose, it suffices to check 
that the eigenvalues corresponding to the reduced density 
matrix are all positive semi-definite. That requirement 
is succinctly contained in the pair of inequalities p\pz ^ 
p 2 _ and pi |p+ 1- In the event, it turns out that the 
latter inequality is all that is needed here, because the 
former condition is already satisfied for all values of the 
parameter B. 

A note of caution with regard to predictions on the 
basis of spin chains with finite length: such finite chains 
are prone to uncharacteristic behaviour, particularly in 
the vicinity of a quantum phase transition. As an exam- 
ple, let us consider how the logarithmic negativity with 
respect to a bridge pair of spins varies as a function of 
the spin chain's length. Such a scenario is depicted in 
Fig. (|13fl , with the parameter B chosen to lie within the 
range of values where the logarithmic negativity is non- 
vanishing in the thermodynamic limit (see Fig. iJHj); 
B = 0.9875 to be precise. 



Fig. (|13fl points to a remarkable finite-size effect: the 
diagram is composed of three separate curves, the middle 
one representing chains of odd lengths and the other two 
together making up the even lengths. Furthermore, one 
of the curves (representing lengths where N + 2 is a mul- 
tiple of four) stays identically zero to start off with, until 
it abruptly rises to converge with the other two curves. 
The diagram sees the logarithmic negativity converge to 
a value around 0.0226, which is in complete agreement 
with the plot of Fig. I|12f> at the same value of B. For 
all intents and purposes, N — 700 can therefore already 
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FIG. 13: The logarithmic negativity with respect to a bridge 
pair of spins is shown as a function of the spin chain's length, 
with the magnetic field strength parameter B set to 0.9875. 
The graph indicates a marked finite-size effect, which van- 
ishes in the thermodynamic limit. The middle curve repre- 
sents chains with an odd number of spins, while the other 
two curves make up the even chain lengths. All three curves 
converge to a common value as the chain length increases. 



be regarded as representing an infinitely long chain for 
this particular value of B. Predictably, as the coupling 
parameter is steadily increased to B c = 1, the three 
curves converge to an increasingly higher value (reach- 
ing 0.5711... at B c — 1), and the rate of convergence 
increases simultaneously. However, when we increase 
the coupling parameter further still, another surprising 
finite-size effect occurs. Precisely where the coupling pa- 
rameter surpasses the critical point B c = 1, the bottom 
curve of Fig. (|13l) (now strictly positive for all N) expe- 
riences a dramatic jump and approaches the convergence 
with the other two curves from above. The size of this 
phenomenon diminishes with increasing chain length and 
vanishes all together for infinitely long chains, as is evi- 
dent from Fig. I|12[l. So it would actually appear that the 
much sought-after critical behaviour in the entanglement 
measure turns out to be a finite-size effect that vanishes 
in the thermodynamic limit! 

In summary, we embarked upon the preceding investi- 
gation in quest of a simple order parameter that would 
allow us to predict the occurrence of an infinite localis- 
able entanglement length. A combination of two-point 
correlation functions motivated by quantum information 
science, the logarithmic negativity was deemed to be a 
suitable candidate for such an order parameter, but our 
findings have since led us to conclude that this measure 
cannot, in general, predict the behaviour of the localis- 
able entanglement length. This further substantiates the 
belief that the localisable entanglement length in transla- 
tion invariant systems is a novel concept that transcends 
mere two-point correlation functions. 



VII. QUANTUM COMPUTATION 

We have already seen how, in an optical lattice, for 
a single triangle, one can manipulate, by varying suit- 
ably the tunneling and/or the collisional couplings, the 
interaction terms of the form crjcrj +1 + 0jf<T* +1 , cr|cr| +1 , 
o-jO-j +1 (J z J+ 2 a*(aj +1 a* +2 +o-y +1 o-j +2 ). These interac- 
tions, up to common single qubit rotations, are equivalent 
to quantum gates where the states 1 1) and | j) represent 
the logical states | 0) and 1 1) of the computation. From 
the above interactions, one can obtain SWAP, controlled- 
Phase (CP), controlled-controlled-Phase (C 2 P) and 
controlled- SWAP (cSWAP) respectively. In the follow- 
ing subsection we shall see how this is achieved. In par- 
ticular, the CP and C 2 P gates are produced by manip- 
ulating one of the two optical lattice modes while the 
other remains with a high amplitude, corresponding to 
zero tunneling. The exchange interaction (SWAP) is pro- 
duced by lowering the barriers that trap both the | 0) and 
1 1) states as seen in <|3.1|) . 

We shall then proceed to show how these gates can 
be combined to perform quantum computation on one- 
or two-dimensional structures, such as those described 
in section without the need for targeting specific lat- 
tice sites i.e. we will only control fields that are applied 
globally, to the whole system [rj Ifjl IBH IH1 l5fj| . 

There are two key ideas in targeting operations to spe- 
cific qubits by global addressing. Firstly, we employ a 
specific qubit as a pointer and, secondly, we consider the 
effect of double wavelength fields for addressing alternate 
triangles. We expand on both of these ideas in subse- 
quent subsections. 

A. Quantum Gates from Interactions 

In order to create quantum gates from the interactions 
that we have, let us define 

Ai= [ \ {l) dt 
Jo 

where = B 13. 1|) . This encapsulates the variation 
with time of the coupling strengths as various optical 
lattices are applied. The results should always be taken 
modulo 27T, due to the exponentiation procedure. The 
time, T, will be the same for each A, for a given gate but 
can, naturally, be different for the different gates that we 
wish to create. 



Gate 


A 


Ax 


A 2 


A 3 


A 4 


CP 


7T 
4 


7T 

4 











C 2 P 


7T 
8 


7T 

8 





7T 

8 





SWAP 








7T 

2 








cSWAP 








7T 

4 





— 7T 

4 



Indeed, the above table shows how to create the quan- 
tum gates from the different terms of the Hamiltonian 
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lpTT)l . As presented, both SWAP and cSWAP generate 
an extra phase gate. This effect can also be negated, but 
has been left out for simplicity. 



B. Organisation of the Computer 

We start with a one dimensional ladder of triangles, as 
shown in Figure [S] and again in Figure HlT a'l. although 
this can easily be extended to two dimensions. We label 
one horizontal line as the register array, and the other 
as the auxiliary array. The register array will contain 
the qubits that we perform the computation on while the 
auxiliary array will facilitate the transport of the pointer. 
The whole system is initialised such that every qubit is 
in the | 0) state. We require a pointer to help us localise 
operations. This is achieved by retaining single qubit 
control over a single lattice site, so we can flip that qubit 
into the 1 1) state. 



(a) a 




(c) 



Auxiliary 

Register 

Auxiliary 

Register 

Auxiliary 

Register 



Pointer 



FIG. 14: (a) The one dimensional chain of triangles divided 
into a row of auxiliary and register qubits. All are initialised 
as a | 0} apart from the pointer, which is in the state | 1). 
The interactions required to transport the pointer are also 
indicated, a and j3 exchange interactions can be controlled 
separately, (b) Interactions required to create the CP gate, 
yielding a z on a single qubit. (c) Interactions required to 
create C 2 P gate on alternate triangles. 



C. Moving the pointer and qubits 

We have specified how certain gates can be created 
when the pointer is next to a specific qubit, or when 
the pointer and two specific qubits are all located at the 
vertices of the same triangle. However, we need to know 
how to move the pointer and the qubits so that they can 
interact as we want them to. To achieve this, we use the 
SWAP operation applied to alternate qubits. 

Moving the pointer is a relatively simple matter if we 
have a system of SWAPs available. We just switch be- 
tween the two SWAP-ing modes along the top of the 
chain, denoted by a and (3 in Figure H~4T a) . This allows 
us to move the pointer to any arbitrary position on the 
chain. Exactly the same idea can be applied to the regis- 
ter qubits, provided they are separated by an even num- 
ber of qubits. If they are separated by an odd number of 
qubits, then applying SWAPs just causes the two qubits 
to move together, with constant separation. To avoid 
this, we have to use the pointer to apply a controlled- 
SWAP to one of the qubits so that it becomes separated 
from the other qubit by an even number. We have al- 
ready shown one method for generating this cSWAP, us- 
ing the natural Hamiltonian of the system. It can also 
be built by standard gates, which we shall see how to 
construct in the following. 



D. Superlattices 

In order to activate certain Hamiltonian terms on al- 
ternate triangles, we need to employ the idea of superlat- 
tices. These are obtained by superposing on the trapping 
potentials standing wave fields with a different period. 
The idea is illustrated in Figure ED To generate these 
superlattices, it is not necessary to have a large set of 
lasers with different wavelengths. Instead of setting them 
up in direct opposition, it is possible to create standing 
waves with varying periodicity by introducing an angle 
between them |57l l58j . This angle determines the period 
of the standing wave, di given by 

A 



2sin(cV2) 



The basic idea is that if we can create a CP gate be- 
tween lattice sites along the lines shown in Figure ITlT bh 
then this acts as a a z gate being applied only to the qubit 
adjacent to the pointer. If we can also create a C 2 P gate 
on alternate triangles, for example around the triangles 
shown in Figure HlT cl . then this is just like creating a 
CP gate between the two qubits on the same triangle as 
the pointer without affecting at all the rest of the regis- 
ter qubits. These actions are the main building blocks 
for performing universal quantum computation with our 
setup. 



The required manipulation of potentials for triangular 
lattices demands the activation of tunnelings along cer- 
tain sites, while it should be prohibited along other ones. 
For example, to create the alternating SWAP that we 
require for moving the pointer (Figure a)), we need 
to ensure that we don't activate couplings between the 
register and the auxiliary arrays. We also only want to 
activate couplings on every other horizontal line. Hence, 
we specify that we require a potential of the form 

V K = cos(kx) sin | — = | sin ( — % — kx \ sin | — = + kx | , 
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FIG. 15: Illustration of how a laser with double the wave- 
length of the trapping potential can activate alternate cou- 
plings, as required for the SWAP procedure (for example). 

taking the origin to be located on one of the register 
qubit lattice sites. This can be expanded as a series of 
sine functions, each of which can be created by a pair of 
lasers resulting in the desired pattern, given in Fig. 1161 
The determination of orientation and angle for the laser 
pairs is demonstrated in @ using the simpler example of 
a square lattice. 




FIG. 16: The potential offset that activates tunneling tran- 
sitions between alternate lattice sites. The darker areas indi- 
cate lowering of the potential barrier, while the lighter areas 
indicate its increase. The straight lines indicate constant po- 
tential. The qubits are located at the vertices of these lines. 



E. Universal Quantum Computation 

So far, we have shown how to create a certain set of 
interactions as we require them - CP, C 2 P, SWAP and 



cSWAP. If we can also apply single qubit rotations to al- 
ternate lines of the lattice, then we can apply sequences 
such as U, CP and U\ which have no effect if the aux- 
iliary qubit is in the | 0) state because U and W cancel, 
but applies Ua z U^ to the qubit targeted by the pointer. 
We are then able to create, by this process, any single 
qubit rotation on a given qubit. Together with the CP 
gate (or VSWAP gate, which is generated using the same 
procedure as for the cSWAP, for only half the time), it 
constitutes a set of universal gates for quantum compu- 
tation. 

The single qubit rotations, created by Raman transi- 
tions on alternate rows of the lattice, can, again, be gen- 
erated by a standing wave of double the period, so that 
the 0's of potential are localised on the rows that we don't 
want to experience the rotation. The lasers that create 
this standing wave need to be of a different wavelength to 
all the others to avoid undesired transitions taking place. 



F. Error Prevention and Correction 

Errors are the bane of any implementation of quan- 
tum computation, and there is much research on min- 
imising their effects with better qu ality ap paratus and 
by using decoherence free subspaces [5|j, |6(| • It is known 
that once we are able to perform operations with an er- 
ror that is below a certain threshold (assuming any errors 
leave the qubits within the computational subspace) , it is 
possible to perform arbitrarily accurate quantum compu- 
tation through the process of concatenation of quantum 
error correcting codes [6l| , leading to fault tolerance [§2 • 
It has recently been shown that the same can be done for 
globally controlled systems such as the one presented here 

iHiia. 

The first thing to note is that the auxiliary array al- 
ways remains in a classical state. This instantly protects 
it against a z errors, since classical states are eigenstates 
of er 2 . We can also protect against a x errors by constantly 
measuring the state of the qubits via the Quantum Zeno 
effect. Note that, if we measure the auxiliary array in 
the 1 1) state, we risk losing the pointer. We choose, in- 
stead, to measure only every other qubit on the auxiliary 
array, not including the pointer. This at least provides 
an indication of whether or not there has been an error 
on the pointer. 

We can also employ this idea to provide measurement 
at the end of the computation. We achieve this by per- 
forming a cSWAP around a triangle so that we move a 
qubit from the register array onto the auxiliary array, 
next to the pointer. We are then able to perform the 
measurement. 

Finally, the qubits on the register array can be en- 
coded into a quantum error correcting code. Sufficient 
parallelism can be generated within the array to allow 
for error correction on every block of encoded qubits (see 
0, |5f| for more details). Instead of correcting the er- 
rors by making a measurement and then performing the 
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relevant correction, these steps are performed by an algo- 
rithm that uses controlled-NOT gates to correct for the 
errors, as determined by the error syndrome stored on 
ancilla qubits. The problem with this method is that to 
perform the syndrome extraction (placing on auxiliary 
qubits information about the errors), the ancilla qubits 
must be in a well-known state (|0)). After the error 
correcting phase, these qubits are left in an unknown 
state. Hence, it is necessary to either reset the state of 
the qubits, or have a large supply of fresh qubits. In the 
optical lattice set-up, this second option is quite sensi- 
ble. We have so far restricted computation to a single 
two-dimensional plane. However, in a three dimensional 
lattice there are many of these planes, all of which contain 
qubits in the | 0) state. If we perform our computation on 
a single plane, then this plane can be moved through all 
the other planes by a series of SWAPs between alternate 
planes. We can therefore access this large supply of fresh 
qubits for the purposes to error correction without the 
need to perform measurements during the computation, 
thus simplifying the experimental implementation. 

VIII. CONCLUSIONS 

In this paper we presented, initially, a variety of differ- 
ent spin interactions that can be generated by a system of 
ultra-cold atoms superposed by optical lattices and initi- 
ated in the Mott insulator phase. In particular, we have 
been interested in the simulation and study of various 
three-spin interactions conveniently obtained in a lattice 
with equilateral triangular structure. The possibility to 
externally control most of the parameters of the effective 
Hamiltonians at will renders our model as a unique lab- 



oratory to study the relationship among exotic systems 
such as chiral spin systems, fractional quantum Hall sys- 
tems or systems that exhibit high-T c superconductivity 
Furthermore, unique properties related with the 
critical behaviour of the chain with three-spin interac- 
tions has been analysed (see also where the two- 
point correlations, used traditionally to describe the crit- 
icality of a chain, seem to fail to identify long quantum 
correlations, suitably expressed by a variety of entangle- 
ment measures 45]. In particular, analysing the logarith- 
mic negativity, indicates a possible connection between 
the localisable entanglement length on the one hand and 
entanglement properties of closely spaced spins on the 
other. In addition, suitable applications have been pre- 
sented within the realm of quantum computation 12?! ] 
where three-qubit gates can be straightforwardly gener- 
ated from the three-spin interactions. 

In conclusion the three-spin interactions generated in 
an optical lattice offer a rich variety of applications in 
quantum information technology as well as in solid state 
physics, worth pursuing further theoretically and exper- 
imentally. 
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